Analysis and Geometric Optimization of Single Electron Transistors for Read-Out in 

Solid-State Quantum Computing. 
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The single electron transistor (SET) offers unparalled opportunities as a nano-scale elec- 
trometer, capable of measuring sub-electron charge variations. SETs have been proposed for 
read-out schema in solid-state quantum computing where quantum information processing 
outcomes depend on the location of a single electron on nearby quantum dots. In this paper 
we investigate various geometries of a SET in order to maximize the device's sensitivity to 
charge transfer between quantum dots. Through the use of finite element modeling we model 
the materials and geometries of an AI/AI2O3 SET measuring the state of quantum dots in the 
Si substrate beneath. The investigation is motivated by the quest to build a scalable quan- 
tum computer, though the methodology used is primarily that of circuit theory. As such 
we provide useful techniques for any electronic device operating at the classical/quantum 
interface. 
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I. INTRODUCTION 

The single electron transistor (SET)^*^ has presented 
itself as a very interesting device in the field of classical 
computing. Systems using SETs as logic devices have al- 
ready been proposed^. SETs also present the possibility 
of utilizing the spin orientation of electrons in the cur- 
rents passing through them. This leads to their use in 
the field of spintronics^, opening a plethora of possibil- 
ities in which this extra degree of freedom can be used to 
encode information. A SET can also be used as a highly 
sensitive electrometer capable of measuring charge varia- 
tions as small as 10~^ge/\/Hz which is close to the quan- 
tum limit^. Our research is motivated by the use of the 
SET to perform solid-state quantum computer read-out, 
though the techniques employed here are applicable to 
any application using a SET as an electrometer. 

The field of quantum computing began in the 1980s 
when researchers began to theoretically investigate how 
to apply quantum mechanics to information processing. 
The first application envisioned for such quantum infor- 
mation processing was as a physics simulator, presented 
by Feynman in 1982 when he noted that classical systems 
can not simulate quantum systems efficiently^. In 1985 
Deutsch developed a quantum algorithm.^ demonstrat- 
ing the potential for quantum information processing in 
its own right. However it was not until the discovery by 
Shor in 1994 of a quantum algorithm for polynomial-time 
factorization of prime numbersSi that the possibilities for 
quantum information processing other than simulating 
physics began to capture the mainstream imagination. 
Shor's work provided much impetus to the field due to the 
ramifications to RSA encryption. Grover developed an- 
other quantum algorithm^ in 1997 for improved database 



searching. This has spurred further interest in the de- 
velopment of architectures for quantum computers, and 
other useful quantum algorithms. Since that time the 
race to build a working quantum computer has taken off 
across the globe. 

The essential difference between classical computers 
and quantum computers is the replacement of the dualis- 
tic on/off bit (0 or 1), with the quantum-bit (qubit) which 
is capable of being in both the on and off state simulta- 
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Qubits may also make use of the purely quantum me- 
chanical property of entanglement, which is the source of 
the unique computational power of quantum information 
processing's. The operation of any quantum computer 
requires that qubits can be individually 'rotated' (i.e. the 
values of a and [3 can be manipulated for each qubit) and 
that two qubit interactions can be controlled (to allow 
for multi-qubit logical operations such as CNOT). While 
these are the essential qubit requirements for universal 
information processing, an actual physical implementa- 
tion requires further considerationsii. 

One of the biggest difficulties in quantum information 
processing is that the encoded information is extremely 
fragile. A delicate balance between optimizing the de- 
sired entanglement between qubits, and minimizing the 
undesired entanglement between the qubit and its envi- 
ronment must be achieved. Any interaction a qubit might 
have with its environment causes the information stored 
on the qubit to "leak" out. This is referred to as decoher- 
ence. Until the development of quantum error correction 
algorithmsiSii^ it was believed by many that this deco- 
herence would prevent the construction of any practical 
quantum computer. Even though constructing a quan- 
tum computer is still at the limits of current fabrication 
technologies, these error correction algorithms show that 
a working quantum computer is (in principle) possible, 
as long as the error rate on the qubits can be kept below 
some threshold (a commonly used estimate^^ii is 10~^). 
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In this paper we are concerned with the read-out stage 
of a siUcon based quantum computers. Read-out con- 
stitutes projecting the quantum information stored in a 
qubit onto the sub-space of classical information of the 
macroscopic measuring device. Qubit read-out will, in 
general, be probabilistic. The read-out process for the 
state a|0)-t-/3|l) will determine a or 1 with respective 
probabilities of jap and |/3p. We consider the read-out 
process for two proposed Si:P solid state quantum com- 
puters (SSQC), the Kane SSQC^^ (FigCJ) and the charge- 
qubit SSQC^^ (Fig|21). Much theoretical work has been 
done on optimizing the read-out process based on how 
to encode the quantum information and which degree of 
freedom of the qubit we choose to measur©i2*iS*i^. Most 
of these proposals essentially require the determination 
of the location of a single electron for read-out. We have 
focused on the SET design in order to optimize it as an 
electrometer for determining this electron location. 

The Kane quantum computer (Fig^) uses the nuclear 
spin orientation (with respect to an external global mag- 
netic field) of a single P dopant in silicon as its qubit. 
The relaxation time of the nuclear magnetic moment of 
P in silicon was extensively studied in bulk resonance ex- 
periments in the 1960sS£, and is known to be in excess 
of 10 hours at low temperatures. This makes the sys- 
tem highly attractive with respect to decoherence issues. 
Since the nuclei we wish to store our information on are 
so well isolated, read-out is extremely difficult. 

This difficulty in read-out can (in principle) be over- 
come using spin to charge transductioni^ as follows. A 
P atom has five electrons in its outermost shell. Four 
of these are used to bond to the surrounding Si lattice, 
leaving a single valence electron loosely bound to the P 
nucleus. The spin orientation of this electron can be cor- 
related with the nuclear spin through the hyperfine inter- 
action. In this way the electron can be used to address 
and infer information about the nuclear spin state. The 
strength of the hyperfine interaction is determined by the 
electron's wavefunction at the nucleus. By carefully con- 
trolled voltages on gates placed ontop of the Si substrate 
(A gates in Fig^l, the electron/nucleus overlap can be 
controlled to bring the P nuclear magnetic moment into 
and out-of resonance with an external, global, oscillating 
magnetic field. This allows for single qubit rotations. 

The electron's wavefunction extends across many lat- 
tice sites. By placing appropriate voltages on surface 
gates placed between qubits (J gates in Fig^, the wave- 
function overlap of neighboring P dopants' electrons can, 
in principle, be controlled, allowing for controlled two 
qubit interactions. In this way multi-qubit logic gates 
can be achieved. In practice such control will require 
precise atomic placement of the P donorsSLSS. This use 
of the valence electron as a 'messenger' between qubits 
can only occur if the decoherence time of the electron 
spin system is long enough to allow for information pro- 
cessing operations to occur. At sufficiently low temper- 
atures, electron dephasing times of greater than 60ms 
have been reported^^, which is more than ample. The 




FIG. 1: Schematic of the Kane quantum computer. The 
qubit states are encoded in the orientation of ^'^P nuclear 
spin. Quantum information is mediated by the donor valence 
electrons. A gates control single qubit rotations via the hy- 
perfine coupling between a valence electron and its nucleus. 
J gates control qubit interaction. Read-out is via a SET (not 
shown) in the vicinity of a P donor in a known state. 



Kane SSQC read-out (Fig|21(a)) and allied electron-spin 
proposalsi^iiSiiS consists in transferring the valence elec- 
tron of a P dopant to a neighboring P-dopant (forming a 
doubly occupied D" state on the latter site). Originally 
this scheme involved transfer by a DC field. To avoid 
ionization of the D~ state a resonant transfer technique 
has been proposed24. This D~ state can occur only if the 
two P donors have opposite nuclear spins. This is due to 
carefully correlating the valence electron's spin to the nu- 
clear spin by the hyperfine interaction, and the D~ state 
being Pauli limited. Spin to charge transduction is thus 
achieved, converting the difficult problem of single spin 
detection to charge detection. In this way, by preparing 
a read-out qubit in a known spin orientation, the ori- 
entation of the qubit being measured can be inferred. 

The charge qubit SSQC (Fig[21) also uses P dopants 
in Si for qubits, however a P-P"*" pair constitutes a sin- 
gle qubit. In this way a single valence electron is shared 
between the dopants. The |0) and |1) states are now en- 
coded onto the location of the shared electron. Read-out 
consists of determining this location (Fig|31b). Geome- 
tries have also been proposed that may be sensitive to 
the superposition basis, by either adding an additional 
ionized dono»2^, or placing the SET centrally between 
the donors2&. 

Single qubit rotations for the charge qubit are achieved 
by varying voltages on control gates above the silicon sub- 
strate, in order to force the shared electron from the left 
(|L)) to the right (|i?)) dopant nuclei. Multi-qubit opera- 
tions can be achieved by carefully placing the P-P+ pairs, 
such that the dynamics of one qubit becomes dependent 
on the other. The silicon band-structure implies atomic 
precision placement of the P-P+ pairs is required^^. 

The disadvantage of the charge-qubit is that it is 
more susceptible to decoherence by local charge fluctua- 




FIG. 2: Schematic of the charge-qubit quantum computer. 
The qubit states are encoded in the location of a shared P- 
donor valence electron. S gates control the symmetry of po- 
tential wells about donors. The B gate creates a barrier 
potential to control tunneling between the donor pair. Read- 
out is via a SET (not shown) used to determine the location 
of the shared electron on the P-P"*" donor pair. 



tions due to the long range of the electromagnetic force. 
This has been the subject of a number of theoretical 
investigations^SiSS. Current estimates put the decoher- 
ence rate in the order of lOnsSS. This short decoherence 
time will still allow for gate operations to occur due to 
the faster operation time (estimated to be of order 10 
P3i&). 

From the above description of both SSQC types it is 
apparent that the final read-out signal for both SSQCs is 
dependent on the location of a single electron on subsur- 
face P donors. Hence by using SETs to monitor the local 
electrostatic environment we can perform SSQC read- 
out. Though our investigation is motivated by this spe- 
cific application of SETs, the analysis technique we follow 
is general enough to be useful to any field intending to 
use the SET as an electrometer. 

Our research is motivated by experiments in which sin- 
gle charge transfer between clusters of P-donors was mon- 
itored by two SETs^SiSi (Fig|3I). We use the structures 
of the experimental SETs as a starting point for our in- 
vestigation into SET geometry optimization. In keeping 
with the experimental architecture, our simulations con- 




(=1 Kane qjbll f'» cnai^e qubll 



FIG. 3: Initial and final states of the read-out process for the 
Kane and charge-qubit SSQCs. 



3 




FIG. 4: SEM image (courtesy of F. Hudson et al, UNSW) 
of a typical experimental device providing motivation for the 
SET geometries under consideration. 

sists of donors placed 20nm below a 5nm silicon oxide 
layer. We also consider clusters of donors initially (see 
section llV A|l . rather than single dopants, with the single 
dopant case considered in section IVDI 

We begin our discussion by presenting some introduc- 
tory SET theory, along with a description of the induced 
island charge as a measure of SET sensitivity. The inter- 
play between quantized electrons, and continuous electric 
fields in relation to understanding mesoscopic circuitry is 
then discussed. The initial and final states of the read- 
out process in both types of SSQC is then described in 
detail. In this paper we extend our previous descrip- 
tion of SET sensitivitji^ to calculate both the charge 
induced on a SET island, and the current variations that 
can be expected through the device due to charge trans- 
fer on nearby quantum dots (QD) (i.e. the nuclear or 
charge qubits). We then present a technique for calcu- 
lating these currents in the steady state, based on a rate 
equation and geometric considerations of the SET. The 
experimental setup which motivated our research, and 
the techniques we use to determine the capacitance ma- 
trix of our system are then presented. Finally we present 
the results of simulations in which we varied the SET 
geometry in order to optimize its sensitivity as an elec- 
trometer. 



II. SINGLE ELECTRON TRANSISTORS 

As electronic devices approach the nanometer scale, 
the sea of continuous electronic energy levels normally 
present in metals becomes discrete and parts of the device 
begin to behave as low dimensional quantum systems. 
When this occurs, the passage of single electrons through 
a tunnel barrier can be controlled. By joining two tunnel 
barriers in series (effectively creating an isolated island) , 
it is possible to create a transistor, the current though 
which is determined by single electron tunneling events. 
The controlled tunneling of single electrons comes about 
due to the Coulomb blockade present on the SET island 
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FIG. 5: Discrete energy levels on a SET island. For sin- 
gle electron tunneling to be observable, the energy spacing 
between levels must be greater than thermal fluctuations 
{AE > ksT). (a) First unoccupied energy level on the is- 
land is greater than the Fermi energy (£"/) of the source pre- 
venting tunneling through the potential barrier U. (b) The 
energy levels of the island are shifted by the application of 
a nearby potential, allowing an electron to tunnel from the 
source. The field effect of the excess electron on the dot pro- 
duces a Coulomb blockade. This prevents further electrons 
from tunneling onto to the island. 



(FiglSl) . Any electron tunneling onto the island faces find- 
ing a particular energy level to inhabit, rather than the 
conduction band continuum normally occurring in met- 
als. Tunneling phenomena of single electrons may be 
observed and controlled by varying the height of the is- 
land's energy levels through the application of a nearby 
electric potential (e.g. the gate of the transistor). Ther- 
mal fluctuations in the device must be small enough not 
to override the resolution of the discrete energy levels, 
i.e. /S.E > ksT, where AE is the energy level spacing of 
the island. 

The Coulombic repulsion of electrons present on the is- 
land prevents any further electrons from tunneling onto 
it. Though the field produced by a single electron is 
quite small, the effect of this field on the device is in- 
versely proportional to the size of the island. Hence in 
the nanometer regime the field effects can become quite 
large. This prevention of further electrons moving onto 
the island is referred to as the Coulomb blockade. Only 
when an electron tunneling to the island has an energy 
equal to or greater than one of the unfilled discrete energy 
levels of the island is it able to overcome the Coulomb 
blockade. By tuning the energy level on the island via 



the transistor gate, the energies for having n or n -I- 1 
excess electrons on the island become degenerate. Elec- 
trons will now able to tunnel sequentially from source to 
drain, creating a measurable current through the SET. 

The behavior of the SET can be characterized by 
sweeping the gate voltage for a variety of source-drain 
(transport) voltages. Mapping the current through the 
SET for each point in this gate-transport voltage plane 
produces regions of zero current in which the tunneling 
is suppressed, and the current is essentially zero. These 
regions correspond to Fig|Sl (a), and have a characteris- 
tic diamond shape. Each successive Coulomb diamond 
(in the gate voltage direction) represents a region of the 
parameter-space in which the charging energy of the sys- 
tem is minimized by adding a single electron to the island. 
By increasing the transport voltage, higher charge states 
become available in the tunneling process, reducing the 
efficacy of the Coulomb blockade, as shown in Fig|Sl 




FIG. 6: Normalized current mapped across the gate-transport 
(source-drain) voltage plane (also normalized). The charac- 
teristic Coulomb diamonds show stable regions of zero cur- 
rent. Varying the gate voltage moves the SET island through 
successive states of excess electrons, while increasing the 
transport voltage reduces the Coulomb blockade, by intro- 
ducing higher order charge states. 



A. The SET as an Electrometer 

As described above, the variation of the potential ap- 
plied to the SET gate varies the energy level on the is- 
land, changing the strength of the Coulomb blockade. In 
an analogous fashion, any charge on a nearby QD will 
vary the island's potential. We can associate a capaci- 
tive coupling between the island and the nearby QD. The 
potential variation on the island is a continuous variable 
based on the island-QD capacitive coupling. From the 
definition of capacitance (C = Q/V) we can associate 
this potential variation with an induced charge Sq, which 
may be a fraction of an electron charge. 

The greater this Sq, the stronger the variation in the 
energy levels on the island. As such, this Sq is commonly 
used to describe the sensitivity of the SET to the system 
being measured. The actual signal measured however is 
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T^ + T^ + tF + T^ -9e. (2) 



Inspection of Equation Q tells us immediately that the 
charge induced on the island due to the QD is simply 



(3) 



FIG. 7: Circuit diagram of the SET-QD system. Isolated 
regions with discrete energy levels are represented by solid 
boxes. The dotted outline indicates the region considered as 
the SET island. 



a current through the SET which we calculate in section 

cni 



B. The Induced Island Charge 

We present two methods of determining the charge in- 
duced on an island for the simplest system of a SET cou- 
pled to a single QD. We extend the technique to a two 
QD system for our analysis, though present the single 
QD case here for clarity. The motivation for developing 
two approaches is that the first is a common method in 
the literature, and is very simple, while the second al- 
lows us to determine the energy variations of the entire 
system due to tunneling events. These energy variations 
can then be used to determine current variations in the 
SET. Comparison of the two techniques indicated that 
both gave equivalent results for our simulations. 

The circuit diagram for the system being considered is 
shown in FigEI 



1. Ratio of QD's Total Induced Charge 

Suppose we have a QD charged with n excess electrons 
(g = —nqe) coupled to a SET. Let the self capacitance 
of the QD be Cix;. This self capacitance is simply the 
summation of the QD's capacitive coupling to each of 
the objects in its environment. In our model we assume 
the capacitive coupling to defects in the lattice (dangling 
bonds, local impurities etc.) is negligible so that we have, 

Cis — Cu + Cis + Cid + Cig , (1) 

where the roman subscripts refer to the SET's island, 
source, drain and gate. Conservation of charge necessi- 
tates that the charge on the QD must induce a charge of 
nqe in its environment. Hence we can write the manner 
in which this charge is shared across the environment as. 



2. The Electrostatic Energy Variation of the System 

This second method for determining Jq, though more 
complicated than the previous, has the advantage of al- 
lowing for calculation of energy variations of the system 
due to tunneling events, which is essential for determin- 
ing the current through a SET (see section ITTl|l . 

The electrostatic energy of a capacitor is given hy E — 
q^ /{2C) where q is the magnitude of the charge on one 
of the plates of the capacitor. The electrostatic energy 
of the SET-QD system under investigation is simply the 
matrix equation extension of this, accounting for all the 
capacitances in the system. Any electrodes connected by 
wires to regions outside the system have their potentials 
controlled externally, and so we need only consider the 
charges on any floating electrodes (in this case the island 
and a single QD). The charging energy for the system is 
given by 

E - \q^C^'Q . (4) 

Q is a vector containing the charge on the island and 
the QD and we refer to Ce as the 'energy' capacitance 
matrix, defined below. The charge on the island and QD 
is a combination of actual electrons, and induced charge 
due to voltage differences, 

Q = Q + nqe , (5) 

where n represents the number of electrons on each low 
dimensional quantum system (in this case the island and 
the subsurface QD), while Q represents the continuous 
charge induced on these objects. In matrix form we have. 




Considering only the island for a moment, we can write 
the induced charge from the source, drain and gate elec- 
trodes as 

= CsVs + CdVd + C,gVg , (7) 

where we are referencing all voltages to the island {Vs = 
—Vd assuming the tunnel barriers in our SET are identi- 
cal). For our system then, we can write Q = CcV. We 
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FIG. 8: Plot of difference between two lowest energy levels 
on the island. The blue line is for the quantum dot being 
measured neutral, the red line for a charge of —q^ on the 
quantum dot. Peak to peak distance corresponds to a charge 
of qe, enabling 5q to be determined by the shift between the 
red and blue plots. 



To determine the shift in the system's electrostatic en- 
ergy profile we first determine the two lowest energy lev- 
els of the system for the case of the QD being neutral. We 
set the source and drain voltages to +Vsd/'^ and — V^s(i/2 
(in our simulations Vsd = 0.1 mV) and step through a 
voltage range on the bias gate ±.Vgraax- For each applied 
bias voltage {Vg) we cycle through ±n excess electrons on 
the island and use Equation Q to determine the electro- 
static energy of the system. We then keep the difference 
between the two lowest energies and plot this energy dif- 
ference with respect to the bias gate's voltage (see FiglSl. 

The zeros of this plot are voltages at which the island 
having n, or n -I- 1 excess electrons is degenerate with 
respect to the energy of the system, and correspond to 
maximum currents through the SET. Similarly, the pe- 
riod of the plot corresponds to adding one electron to the 
island. 

To determine 5q we then perform the entire process 
again, this time placing one electron on the QD. This 
will give a similar energy difference plot, but shifted due 
to the dot's influence on the island. Since we know the 
period of the plot corresponds to an electron being added 
to the island, we can thus deduce what fraction of an 
electron corresponds to the shift, giving 5q. 



refer to Cc as the 'correlation' capacitance matrix, 

where the subscript 1 again refers to the QD. This can 
easily be extended to an arbitrary number of QDs. We 
can now write the total charge on the QD and the is- 
land as a function of the electrode voltages and elements 
of the system's capacitance matrix, with the number of 
electrons on each as a free parameter we will vary in order 
to find the minimum charging energy of the system. 

The Ce matrix describes the cross capacitances of the 
QD and the island and is given by, 

where Cis is defined in Equation The island's self 
capacitance (Cis) has an equivalent form. 

While not necessary to the theory, it is worth noting 
that this matrix must be symmetric, which can reduce 
computation time for more complicated systems. We are 
now in a position to calculate the electrostatic energy of 
the system for arbitrary charge configurations on the QD 
and island by simply varying the integers constituting n. 
From this the 6q induced on the island by the presence 
of a charge on the QD will be determined. We do this for 
the case of a single electron being added to the QD, and 
determine 5q by the shift in the system's electrostatic 
energy profile. 



III. CURRENT THROUGH A SET 

Though we have focused on the induced island charge 
{5q) as characterizing the sensitivity of our SET when 
measuring our QD system, the actual signal detected in 
experiments is a variation of current. Having demon- 
strated how to determine the charging energy of our sys- 
tem in section 111 B 21 we now show how to extend the 
method to determine actual currents through the device. 
Our technique closely follows Grabert and Devoret'^^ and 
gives currents in the nA regime, which is in agreement 
with the measured results for SETs in generaU*^. 

To determine the current through the SET we need 
to consider the tunneling rates on and off the island. 
The rate is determined by the energy difference of the 
system between different configurations of electrons on 
the island. We make the assumptions that co-tunneling 
events can be ignored, and assume that only the island 
has quantised energy levels. 

The rate at which electrons tunnel from the electrodes 
to the island is dependent on the number of excess elec- 
trons on the island. Given the island is in a state with n 
excess electrons, and at a temperature T we denote the 
rate at which electrons tunnel from the island (?) to the 
electrode x (either the source or drain) as, 

r". = -1 ^ (10) 

q.Rt exp(^)-l 

We used 4K for our simulations to avoid low temper- 
ature convergence issues in determining the capacitive 
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coupling between the objects. In practice the opera- 
tional temperature of a SSQC would be of order mK, 
though this won't effect the qualitative results of our 
analysis. The parameters of interest are the change in 
energy (Ai?"j) and the barrier tunnel resistance {Rt)- 

We determine AiJ^j by using Equation ^ for the 
charging energy of the system, and considering the work 
done (either on or by the system) in pushing an elec- 
tron through the tunnel barrier. The energy variation 
for Equation (|1U|I is given by: 

AE^^ = E{n - 1) - E{n) + V^q, , (11) 

For the tunnel resitance we use the expression^ 

where W is the width of the tunnel barrier and A is its 
surface area. The variable is the height of the potential 
barrier (taken to be a typical value of 2eV'^^ and Ep is 
the Fermi energy of the island (taken to be 11.65eV at 
4KS). The effective mass of an electron in aluminium 
oxide is taken to be m* = 0.35m(P^. 

To obtain the current through the SET we must con- 
sider the contribution from all possible processes. Since 
we are ignoring co-tunneling processes we need only con- 
sider events that change the number of electrons on the 
island by 1. The current is then given by 

oo 

^ = ^^Y. E VniVl^-^-t^) . (13) 
n— — oo ^—s,d 

where p„ is the probability that the island is in a state 
with n excess electrons. To determine these probabilities 
we consider the master equation 

p„ = (r--i-Hr^i)p„_i (14) 
-(rr, + r:, + rs, + rr>„ 



Equation H14() simply states that the rate of change of 
occupation number must simply be the rate at which 
electrons are entering the state, less the rate at which 
they are leaving. 

In order to solve the problem numerically we must 
truncate the summation. We let the island charge config- 
uration range over ±iV (with N = 10), and Equation H14|) 
becomes a system of 2iV-|- 1 equations. We found = 10 
to be sufficient for convergence, with = 5, = 10 and 
Af = 20 all giving the same results for Equation 1)131) . We 
write these 2A^ -I- 1 coupled equations in a single matrix 
equation. For clarity we make the substitutions 

-pn— 1 . -pn— 1 \ 

^ is ^ ^ id — ^n, 



rrs + C, + rs. + rr, = s„, (i5) 



^ si + ^ di — '-'n • 

Notice that A^n contains the probability for the — A^ — 1 
state, and Cn for the A^ -I- 1 state. We must therefore 
make the further approximation that A-^ = and Cn = 
0. The problem can be simplified further by recognizing 
that we are after the current in the steady state, and the 
Pn are probabilities, hence 

P„ = , (16) 

and 

n 

The matrix equation can therefore be written as (replac- 
ing the final equation with the normalization condition 
of Equation ifTTIl ) 











B 



c 



A — R C 



\ f P-N \ 

J3-« + i 





V 1 



1 1 / V p„ / 



(18) 



r 



Equation (|18|l is a simple matrix equation of the form (|13|l we thus obtain the current through the SET. FigO 
X=YP, which we solve for P, giving all the occupation shows an example of the currents calculated using this 
probabilities. Inserting these probabilities into Equation method. It clearly displays the expected Coulomb block- 
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FIG. 9: Example of currents calculated using the above 
method. Blue plot is for measured QD being neutral, red 
is for the QD having —qe excess charge. The charge on the 
QD induces a shift in the current {51), due to the charge in- 
duced on the island by the QD {5q) behaving in an analagous 
fashion to a voltage variation on the gate. 

ade regions. The two plots arc the currents for the QD 
neutral, and having charge oi —q^. The shift in the plots 
is due to the charge induced on the island due to the QD 
charge. When performing measurements of a SSQC, the 
signal will be the variation in the current, due to the shift 
brought about by the charge on the quantum dot. Our 
results below are thus displayed as the difference between 
the neutral and charged plots. We can maximize our sig- 
nal by making measurements when the SET current is at 
a maximium rate of change with respect to voltage vari- 
ation on the gate (i.e. at maximum transconductance) . 



IV. DETERMINING CAPACITANCES 

A. Experimental Motivation 

We base our SET geometry investigation around recent 
experiments^. The experiment involves two Al SETs 
(with oxide tunnel barriers) atop a 5nm Si02 layer de- 
posited on high resistivity Si. The SETs are coupled 
to a double QD system embedded in the Si substrate. 
Control and barrier gates are used to mediate the elec- 
tron transfer process between the QDs. The gates of 
the SETs compensate the effect of the control and bar- 
rier electrodes on the SETs' islands (See Fig^j). This 
ensures that only uncompensated effects due to electron 
motion between the QDs being measured give a varia- 
tion in the currents through the SETs. Once the full 
capacitance matrix of the system is determined, we use 




FIG. 10: Geometry of experimental motivation for SET ge- 
ometries under investigation. Control and barrier electrodes 
induce charge motion between the clusters. SET gates com- 
pensate the effect of the control and barrier electrodes on the 
SETs' currents, allowing the detection of charge motion be- 
tween the clusters. 

this to calculate read-out signals using the above argu- 
ments. To determine the capacitance matrix of such a 
complicated geometry, it is necessary to use the numeri- 
cal technique of finite element modeling (FEM). In FEM, 
the physics of a system being investigated is calculated 
on a mesh which discretises the geometry of the situation. 
To reduce computing requirements careful consideration 
needs to be given the mesh construction. Finer meshes 
are required in regions where the electric field is expected 
to vary rapidly. Fig llll shows an example of this for the 
single SET single QD system discussed in section ITl B 21 
We use ISE-MESlSSi and in-house computational geome- 
try algorithms to generate the mesh automatically based 
on our input geometery. The experiments motivating our 
research employed twin clusters of 600P donors as a test 
for single electron transfer. We thus modeled our QDs as 
metalic cubes (60 x 60 x 60nm'^) with a surface area of 
approximately the same size as the spheroids determined 
by SRIM^ calculations of the ion implantation process^ 
used to make the clusters. In section IV Dl we investigate 
the change in the measured signal through the SET due 
to measuring single donors rather than clusters. 

B. Capacitance Calculation 

Once the mesh is determined, capacitances are deter- 
mined by performing an AC analysis through the device 
using ISE-DESSIS'^^. This constitutes placing a known 
alternating current between the objects we wish to know 
the capacitance of and solving 

I^A-V + j-oj-C-V , (19) 

/ is the vector containing the current at each of the nodes 
(determined by applying Kirchhoff's current law), V is 
the vector containing the known potential at each node, 
A is the admittance matrix containing the resistance be- 
tween each node, w is the frequency of the AC. We use 
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FIG. 11: Example of mesh variation based on geometry. The 
colours represent the electrostatic potential due to an excess 
electron placed on a subsurface quantum dot. Finer meshes 
near metallic objects are generated automatically to capture 
the rapidly spatially varying electric field. 

a frequency of IMHz, but the model is not sensitive to 
this parameter since we have a pure Si substrate. Finally 
C is the capacitance matrix containing the capacitances 
in the device. It is worth noting that some properties 
of the capacitance matrix can serve as a check for the 
correct behavior of the simulation. Firstly the elements 
of the capacitance matrix (Cq/j) denote the capacitance 
between the ath and /3th objects hence Cap = C/3cc Di- 
agonal terms (Cqq) denote the total self capacitance of 
the ath object. Since the charge on the ath object in- 
duces a charge of opposite sign on all other objects: 

— ^ Cap = Caa ■ (20) 



V. GEOMETRIES INVESTIGATED 

The effects of misalignment of the QDs to the twinSET 
architecture has been considered in previous studies'*^ us- 
ing the FASTCAP package, though the geometry of the 
SETs was not investigated. Due to the computational 
requirements of FEM, we were not able to perform simu- 
lations of the entire twinSET device that motivated our 
investigations (see section llV A|) . We instead chose to 
model a single SET based on the experimental set-up, 
coupled to two QDs. While this does not change the 
qualitative analysis of our results, we do expect the sig- 
nals for the full twinSET system to be smaller due to the 
larger number of metallic elements the charge we wish to 
measure can couple to. When modeling the Kane qubit 
the SET is coupled to two QDs. Initially the current are 
determined with both QDs neutral. The final configura- 
tion is then with one QD having and the other having 
—Qe- For a charge qubit the SET is also coupled to a two 
QDs, one with the charge of a hole, which changes its 
location (see Figl^)). 

For both qubit types we investigated four aspects of 
the SET geometry. Firstly, we grew the antenna from 




FIG. 12: Geometry of tunnel barrier. The vertical and hori- 
zontal areas contribution to the SET current calculated sep- 
arately and added. 



the island towards the QD to confirm that this improved 
the signal for the experimental device. We then repeated 
this setup with the gate placed between the island and 
the drain, in order to shield the island's antenna from 
the gate. An investigation into the signal variation for 
increasing the size of the island was then performed, as 
this will vary its capacitive coupling to the rest of the en- 
vironment. Finally we investigated varying the overlap 
of the source and drain with the island as this varies the 
tunnel barrier area, implying an increase in the detected 
current variation. The geometry variation can be seen 
in Fig ll3l Fig ll4l and Fig^l For the geometries investi- 
gated the barrier wraps over two faces of the island. We 
treat this by calculating the resistance, width and area 
of the vertical and horizontal parts individually, adding 
the results for the current (see Fig ll2|l . 



A. Results of Growing Antenna 

Fig ll^f a') shows the series of geometries investigated 
for growing an antenna from the island to the two QDs 
representing the measured qubit. Increments of the an- 
tenna length for each simulation was lOnm. The proce- 
dure was carried out for the gate placed beside the drain, 
and for the gate placed between the source and drain. 

Fig ll^f b) and (c) show the variation in current through 
the SET with respect to the gate voltage. The variation 
is for the initial and final qubit read-out states for the 
charge and Kane qubits respectively, for each antenna 
length shown in Fig ll3f a). Note that for the purpose 
of comparison we have kept the spacing between QDs 
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constant for all simulations, although donor spacing for 
the Kane and charge SSQC will most likely be different. 
The results clearly show an improvement in the detected 
signal for the read-out event as the antenna is extended 
towards the system being measured. The signals for both 
the charge and Kane qubit read-out events are of the 
same size, with the different qubit types simply giving 
a shift in the required gate voltage for maximizing the 
detected current variation. This trend was consistent for 
all the simulations performed and henceforth only results 
for the charge qubit will be discussed. 

FigEI(d) shows the effect of placing the gate between 
the source and drain. The size of the current variation is 
the same for both cases (charge qubit data displayed in 
figure) , but the periodicity is greatly increased for the be- 
tween case (note the change in the voltage scale between 
FigEl(c) and FigE](d)). This is due to the self capac- 
itance of the island increasing when the gate is placed 
between the source and drain. Though the signal size is 
not affected, the increase in periodicity is not desirable as 
it would make the current in the SET more susceptible 
to voltage fluctuations in the gate. 



B. Results of Increasing Island Size 

FigEJa) shows the geometries used in a series of sim- 
ulations investigating the effect of increasing the size of 
the island. The width of the island was increased by 5nm 
for each simulation. Though the width of the island was 
more than doubled, no significant change in the current 
variation of the SET due to a read-out event was shown 
(FigEIb)). A slight increase in the periodicity of the sig- 
nal did occur. Of course the island must be small enough 
so that the energy level separation due to charging events 
are greater than the average thermal energy of the elec- 
trons in the system. Our results show that this occurs 
for a wide range of island sizes. As such the size of the 
island is not a significant parameter for consideration in 
SET design. It should be noted however that our sim- 
ulations did not include any stray capacitances to other 
objects in the environment near the SET. In reality we 
expect random impurities in the silicon and silicon oxide 
regions, to which the island will capacitively couple, re- 
ducing the signal size, implying a smaller island size may 
be desirable. 



C. Results of Varying Tunnel Barrier 

Fig El (a) shows the geometry for a series of simula- 
tions investigating the effect of increasing the overlap of 
the source and drain with the island. The overlap was 
increased by 5nm for each simulation. We expect to see 
an increase in the current through the SET as the over- 
lap is increased, as more area is available for electrons 
to tunnel through. This is clearly shown in Fig^] (b) 
which displays the current variation for a charge qubit 





FIG. 13: (a) Geometries for a series of simulations investi- 
gating the effect of an antenna on the SET island. Each in- 
crement of the antenna length is lOnm. Simulations were 
run with the gate beside the drain, and then between the 
source and drain, (b) Current variation for growing antenna 
with SET gate beside drain when measuring a charge qubit 
read-out event, (c) Current variation for growing antenna 
with gate beside drain when measuring a kane qubit read-out 
event, (d) Current variation for growing antenna with gate 
between the source and drain when measuring a charge qubit 
read-out event. 



read-out event for the geometries in Fig^] (a). Al- 
though an increase in current is desirable as it gives a 
larger detectable signal, note that for large overlaps, the 
Coulomb blockade no longer completely suppresses the 
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FIG. 14: (a) Geometries for a series of simulations investi- 
gating the effect of the SET island size. Each simulation 
increased the width of the island by 5nm. (b) Current varia- 
tion for increasing island size when measuring a charge qubit 
read- out event. 

current through the SET. This is shown in Fig ll6l where 
the larger tunnel area simulation always has some cur- 
rent passing though it. This would mean that during 
all qubit operations the SET would have a non-zero cur- 
rent. This would be detrimental to the operation of the 
SSQC, as the back-action of the SET current on infor- 
mation processing qubits would contribute to the overall 
decoherence. We found that simulation 'i' in Fig 1 151 was 
the maximum overlap ( 30nm x lOnm -|-10nm x lOnm) 
that could be used before this occurred. 



D. Results of using Single Donors instead of 
Clusters 

The previous analysis has been based on measuring 
QDs based on clusters of ~ 600 P donors. Based on 
data from SRIM calculations the QDs were modeled as 
60x60x60 nm^ metallic cubes. In order to model the 
single donor case the previous simulations were repeated 
with the original QDs replaced with 6x6x6 nm'^ QDs. 
6nm was chosen based on the Bohr radius of a P donor 
valence electron in Si being ~ 3nm. 

The results of this are displayed in Fig El and show 
that the reduction in QD size has not changed the size 
of the detected signal, merely the required voltage for 



maximizing the signal. This initially seems contradictory 
to the data in Lee et al4i, which found a variation in 
the induced charge for larger QD. This work differs in 
that the variation in the current due to a transfer event, 
rather than simply the charge induced on an island due 
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FIG. 15; (a) Geometries for investigating the effect of in- 
creasing the overlap of the source and drain with the island. 
The overlap was increased by 5nm for each simulation, (b) 
Gurrent variation for increasing overlap when measuring a 
charge-qubit read-out event. 




FIG. 16: Currents through SET for geometries a,i and j in 
Fig llSI Increasing the barrier area suppresses the Coulomb 
blockade. 
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Current Variation (pA) 

FIG. 17: Current variation for charge qubit read-out event 
with SET gate beside SET drain (as in geometry displayed in 
Fig |13l (a')'). The blue lines are the same as Fig |13l (b') with the 
coloured lines for the same simulation with the quantum dots 
being measured being replaced with 6x6x6 nm^ structures 
instead of 60x60x60 nm"^ structures. The plots indicate that 
the signal will be of the same size when performing measure- 
ments with single donors rather than clusters of ~ 600. 

to a two QDs is taken as the cause of the variation in 
the detected signal. Both QDs can capacitively couple 
more strongly to the island^, while the charge variation 
on the island due to an electrons motion between them 
remains unchanged. 

VI. CONCLUSION 

We have investigated the effect of SET geometry varia- 
tions on the detected signal for measurement of a double 
quantum dot system. The two QD system represents 
the read-out process for two types of SSQC qubits, the 
charge and spin qubits. The two qubit types are differ- 
entiated by the initial and final charge configurations on 
the QDs (see Fig|2)l. QDs of size 60x60x60 nm'^ were 
used initially, based on experiments using clusters of P 
donors. We found that growing an antenna to couple the 
island to the QDs being measured greatly increased the 
detected signal for qubit read-out events of both types. 
Variation in qubit types caused no change to the size of 
the signal, and only caused a shift in the required gate 
voltage for maximum current variation. This was found 



to be the case for all geometry variations studied. 

The position of the gate was moved from beside the 
island to between the source and drain. This resulted in 
an increase in the periodicity of the current fluctuations 
through the SET. Since this resulted in no increase in the 
actual signal size, we conclude that placing the gate to 
the side would be more desirable, as the SET would be 
more stable with respect to voltage fluctuations on the 
gate. 

The size of the island was then increased to more than 
twice its original width. This also increased the period- 
icity of the current fluctuations, although the effect was 
negligible. Hence, (ignoring stray capacitances to local 
impurities) the size of the island is not a sensitive param- 
eter in this architecture, as long as the energy spacing of 
the levels on the island remain greater than the thermal 
energy of electrons tunneling through the SET. 

We then varied the overlap of the source and drain with 
the island, effectively increasing the area through which 
electrons can tunnel on and off the island. Increasing 
this overlap increases the signal of read-out events, how- 
ever if the overlap becomes to large, the Coulomb block- 
ade no longer completely suppresses currents through the 
SET. To avoid excess decoherence from the SET when 
measurements are not required, the SET will need to be 
switched on and off. As such, the size of the overlap needs 
to be as large as possible while still allowing for the SET 
to be reliably turned off. For widths of lOnm we found 
that the maximum overlap allowable was 30nm. 

Finally we repeated all the simulations with QD of 
6x6x6 nm^ instead of 60x60x60 nm"^ in order to in- 
vestigate the affect of using single donor atoms instead 
of clusters of ~ 600. The size of the signal we expect 
to see was not reduced at all, with the only effect being 
a shift in the required gate voltage for maximum signal 
detection. 

We wish to thank K. Lee (UNSW) for helpful discus- 
sions with respect to experimental SET geometries. We 
also wish to thank F. Hudson and the measurement team 
at UNSW for the SEM image of the twinSET device. 
This work was supported by the Australian Research 
Council, the Australian government, the US National Se- 
curity Agency, the Advanced Research and Development 
Activity and the US Army Research Offlce under con- 
tract number DAAD19-01-1-0653. 



1 T. A. Fuhon and G. J. Dolan, Phys. Rev. Lett. 59, 109 
(1987). 

2 M. H. Devoret and R. J. Schoelkopf, Nature 406, 1039 
(2000). 

^ K. Uchida, IEEE Transactions Electron Devices 50, 1623 
(2003). 

* J. H. Smet, R. A. Deutsch, F. Ertl, W. Wegscheldar, and 

K. von Klitzing, Nature 415, 261 (2002). 
^ G. A. Prinz, Science 282, 1660 (1998). 
® R. Feynman, International Journal of Theoretical Physics 



21, 467 (1982). 

D. Deutsch, Proc. R. Soc. London. A (1985). 
* P. Shor, Proceedings of the 35th Annual Symposium on 

Foundations of Computer Science (1994). 
8 L. K. Grover, Phys. Rev. Lett. 79, 325 (1997). 
" R. Blume-Kohout, C. M. Caves, and D. Deutsch, Found. 

Phys. 32, 1641 (2002). 
^ D. P. Divencenzo, Fortschritte der Physik 48, 700 (2000). 
^ A. M. Steane, Phys. Rev. Lett. 77, 793 (1996). 
^ P. Shor, Proceedings, 37th Annual Symposiumon Founda- 



13 



14 



tions of Computer Science (1996). 
J. Preskill, arXiv:quant-ph/9712048 (1997). 



^5 B. Kane, Nature 393, 133 (1998). 
" L. C. L. HoUenberg, A. S. Dzurak, C. Wellard, A. R. Hamil- 
ton, D. J. Reilly, G. J. Milburn, and R. G. Clark, Phys. 
Rev. B. 69, 113301 (2004). 
" R. Vrijen, E. Yablonovitch, K. Wang, H. W. Jiang, A. Ba- 
landin, V. Roychowdhury, T. Mor, and D. DiVincenzo, 
Phys. Rev. A. 62, 012306 (2000). 
^® C. D. Hill, L. C. L. HoUenberg, A. G. Fowler, C. J. Wellard, 
A. D. Greentree, and H.-S. Goan, arXiv;qunat-pli/0411104 
(2004). 

R. de Soussa, J. D. Delgado, and S. D. Sarma, Phys. Rev. 

A. 70, 052304 (2004). 

2° D. K. Wilson and G. Feher, Phys. Rev. 124, 1068 (1961). 

B. Koiller, X. Hu, and S. D. Sarma, Phys. Rev. Lett. 88, 
027903 (2002). 

^2 C. J. Wellard, L. C. L. HoUenberg, F. ParisoU, L. M. Ket- 
tle, H. Goan, J. A. L. Mcintosh, and D. N. Jamison, Phys. 
Rev. B. 68, 195209 (2003). 

A. N. Tyryshkin, S. A. Lyon, A. V. Astashkin, and A. Rait- 
simring, Phys. Rev. B. 68, 139207 (2003). 
2" L. C. L. HoUenber, C. J. Wellard, C. I. Pakes, and A. G. 
Fowler, Phys. Rev. B. 69, 233301 (2004). 

A. D. Greentree, A. R. Hamilton, and F. Green, Phys. Rev. 

B. 70, 041305 (2004). 

T. M. Stace, S. D. Barett, H.-S. Goan, and G. J. Milburn, 
Phys. Rev. B. p. 205342 (2004). 

X. Hu, B. Koiller, and S. D. Sarma, cond-mat/0412340 
(2004). 

^® S. D. Barrett and G. J. Milburn, Phys. Rev. B. 68, 155307 
(2003). 

L. Fidichikin and A. Fedorov, Phys. Rev. A. 69, 032311 



(2004). 

^° R. G. Clark et al., Philos. Trans. R. Soc. London, Ser. A 
361, 1451 (2003). 

T. M. Buehler et al. (2004 in preparation). 

C. I. Pakes, V. Conrad, J. C. Ang, F. Green, A. S. Dzurak, 
L. C. L. HoUenberg, D. N. Jamieson, and R. G. Clark, 
Nanotechnology 14, 161 (2003). 

H. Grabert and M. H. Devoret, Single charge tunneling : 
Coulomb blockade phenomena in nanostructures. (Plenum 
Press, New York, 1992). 

D. E. Grupp, T. Zhang, G. J. Dolan, and N. S. Winggreen, 
Phs. Rev. Lett 87, 186805 (2001). 

G. Schon, Quantum Transport and Dissipation (Wiley- 
VCH, Weinhein, New York, 1998), chap. Single-Electron 
Tunneling, 1st ed. 

K. Gloos, P. J. Koppinen, and J. P. Pekola, J. Phys: Con- 
dens. Matter 15, 1733 (2003). 

R. Brenner, Thesis, University of New South Wales (2004). 
Y. Xu and W. Ching, Phys. Rev. B. 5, 4461 (1991). 
^® part of the Integrated Systems Engineering Technology 
Computer Aided Design (ISE-TCAD) semiconductor finite 
element modeling suite., URL www.ise.com 
C. Yang, D. N. Jamieson, C. I. Pakes, S. Prawer, A. S. 
Dzurak, F. E. Stanley, P. Spizziri, L. Macks, E. Gauja, 
and R. G. Clark, Japanese Journal of Applied Physics 42, 
4124 (2003). 

K. H. Lee, A. D. Greentree, J. P. Dinale, C. C. Escott, A. S. 
Dzurak, and R. G. Clark, Nanotechnology 5, 74 (2005). 
SRIM: Stopping and Range of Ions in Matter computer 
model based on J.F.Ziegler, J. P. Biersack and U. Littmark 
- The Stopping and Range of Ions in Solids, Pergamon 
Press, New York, (1985). 



